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An elementary introduction to functional methods and the Schwinger-Dyson equations is pre- 
sented. Emphasis is placed on practical topics not normally covered in textbooks, such as a di- 
agrammatic method for generating equations at high order, different forms of Schwinger-Dyson 
equations, renormalisation, and methods for solving Schwinger-Dyson equations. 
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I. INTRODUCTION 



These notes introduce the derivation and uses of Schwinger-Dyson equations (and nonperturbative methods in 
general) to students who have some familiarity with quantum field theory These methods are based on the functional 
approach to quantum field theory and the path integral, and draw heavily on analogies to statistical mechanics. Thus 
exposure to these concepts will be beneficial to the student. 

The path integral for quantum mechanics is reviewed very briefly and then applied to fields. This leads to a short 
introduction to functional calculus, diagrammatics, and the effective action. The power of the functional method 
is illustrated with a simple demonstration of Goldstone's theorem for the 0(2) sigma model. An aside on the 
functional Hamiltonian and the variational principle follows. We then develop the Schwinger-Dyson master equation 
in a convenient form and derive a diagrammatic method to generate Schwinger-Dyson equations. Applications to ip 4 
theory a contact fermion model, and QED follow. Short discussions of renormalisation, approximation schemes, and 
numerical techniques are included. 

Although these notes are introductory, the discussion will occasionally get technical. Skipping such technical 
sections will do no harm to understanding the material. If you are familiar with things like generating functionals 
and 1PI diagrams then you can just skip to section Ivl 



II. NONPERTURBATIVE QUANTUM FIELD THEORY 



One often hears the word 'nonperturbative' in relation to quantum field theory, and especially Quantum Chromo- 
dynamics. But its meaning is sometimes distorted and it is worthwhile to dwell on it for a moment. 

In the context of quantum field theory, a nonperturbative property is one that cannot be obtained in perturbation 
theory. Notice that this does not mean that a coupling constant is large. For example, the series 

F = 1 + x + x 2 + . . . (1) 
does not converge for \x\ > 1, however one can formally sum it to obtain 



F=- (2) 



and then use this as a definition of F for all x. 
But consider the function 



F = c- 1 / X . (3) 
The Taylor series for this function about the origin is 



F = + + 0+ ..., (4) 

which is not terribly informative. Thus if the only tool one had to examine 'F-theory' were perturbation theory about 
x = 0, no progress could be made. We say that F-theory is nonperturbative. 

There is a common impression that being nonperturbative is akin to a death knell for a theory, consigning it to a 
class of theories that can only be examined with numerical methods. This is overly pessimistic, in the case of F-theory 
a sufficiently acute student might examine it in a new way and determine that it obeys the condition 



(5) 
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This is an easy equation to solve and tells us everything about F, up to a constant. 

Perturbation theory is often the only tool a student learns to compute with quantum field theory. But as we will 
see, many field theories are nonperturbative and other methods must be employed. As with F-theory, a declaration 
that a certain problem is nonperturbative should not be considered fatal. Rather it is an indication that creativity is 
required. 



III. SCALAR FIELD THEORY 



Following long tradition, we shall focus on the properties of scalar field theory in most of these notes. More realistic 
theories will be considered in section IV Dl 

Scalar field theory is a relativistic model of self-interacting spin-zero particles. The 'self-interacting' might be a bit 
odd for students accustomed to particles interacting via vector boson exchange; in this case forces are generated by 
scalar boson exchange, and there are no 'matter' particles. But interesting things can still happen: 'phions' could 
bind into multiphion particles, parity symmetry could be broken, particle interactions could be screened to the extent 
that they disappear at moderate distances. 

A free phion is described with the Lagrangian 

L = J d^x (±d^<f - \mlA . (6) 
Self interactions are introduced via an interaction term L = Lq — J V . This is often taken to be 



(7) 

One also sees interactions of the type V = gip 3 but this theory does not have a finite classical vacuum (as can be seen 
by letting ip be a constant and considering <p < 0), so we do not consider it here. Thus our classical action is 



S = I d d x[ \d^d^ - \m 2 p 2 - ^ 4 ) . (8) 



2 r 2 ur 

Some remarks on the notation: 

(i) We work in d spacetime dimensions. Since (with h = 1) the action is unitless, one obtains 



[mo] = 1 
[<p] = -l + d/2 

[A ] = 4-d (9) 

where units are measured in terms of energy. Thus [mo] = 1 means that this parameter can be specified in GeV. 

(ii) At this stage, the meaning of the parameters mo and Ao is not known. One must compute observable properties, 
such as the mass of a phion molecule or a phion-phion scattering cross section, to enable an interpretation of the 
parameters of the theory. A noninteracting phion obeys the dispersion relation E 2 = p 2 + m$ so it is tempting to 
interpret mo as the phion mass. But it must be remembered that this is only an approximation, interactions will 
change this relationship, and hence the interpretation of mo. For this reason, mo and Ao are called the bare parameters 
of the theory. The physical phion mass, m, is related to mo and Ao in a complicated way. (A similar statement holds 
for the coupling A). 

A famous property of quantum field theory is that calculations generate infinite expressions. The process of 
renormalisation absorbs these infinities in the bare parameters, leaving well-defined finite physical parameters, m 
and A. Renormalisation in the context of Schwinger-Dyson equations and the functional method will be discussed in 
Section [V El For now we note that renormalisation would be required even for a theory with no infinities because, as 
we have stressed, m and mo are not related in a simple way. 

The classical field equations for <p 4 theory are 



(0- 



+ m 2 



(10) 
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This is a nonlinear relativistic field equation that must already be solved with numerical methods or inspired guess- 
work. Equations of this sort can often yield surprising phenomena such as dispersionless waves (solitons). Rich 
nonperturbative phenomena in the quantum realm should be expected! 



IV. FUNCTIONAL METHODS 



A. Introduction 



Following Dirac, canonical quantisation of </? 4 theory proceeds by replacing oscillator strengths in the Fourier 
transform of ip with ladder operators. We summarise this procedure with the expression 

<p-Kp. (11) 
Scattering amplitudes can be computed with the aid the Gell-Mann-Low equation: 

(outlSlin) = m, (12) 

Evaluating the right hand side requires computing matrix elements of products of fields, for example 
(out | / d d x<p 4: (x)\m). This laborious exercise in ladder algebra is greatly simplified with the aid of Wick's theorem: 



T[0(x)0(y)} = N[0(x)0(y)} + (O\0(x)<f(y)\O). (13) 

Besides saving labour, Wick's theorem permits computations solely in terms of field operators, rather than Fourier 
components of fields. Furthermore, practice with Eq. [13]reveals that the algebra of Wick's theorem is the same as that 
of derivatives. These seemingly trivial observations carry great weight because they hint that it is possible to quantise 
field theory solely in terms of the functional properties of fields. This forms the basis of functional quantisation. 



B. Path Integral — Quantum Mechanics 

Functional quantisation takes the Feynman path integral as its starting point. I assume that the reader has been 
exposed to the path integral before, but will nevertheless take this opportunity to review some, possibly lesser known, 
aspects of this method. 

The path integral is interpreted as a sum of phases associated with all possible paths that a system can take from 
initial to final state. In particular, a transition amplitude is written as 

(x f ;tf\xi;U) = J Vxe^ti dtC{± ' x) . (14) 
However, as it stands this expression is meaningless. In fact a better way to write it is 



{x r M^\U) = Y. eiSol{path) ( 15 ) 

where the sum extends over all paths connecting (xi,ti) with (xf,tf). This is traditionally represented by discretising 
the time axis into slices of size e and integrating over all x n — x(ne). Notice however that this means that we are no 
longer summing over all paths but only those which obey the order 1/2 Lipschitz condition: 



\x(t')-x(t)\<k\t' -t^ 2 (16) 

where k is some constant. 

Thus we rewrite Eq. [14] in a more precise form as 
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, N , 

(x / ;t / K;t i > = hm o /n^e*^o^(^x 3 . +1 ) (17) 

J 3=1 

where S c i(xj — > Xj+x) is the classical action evaluated between x(Je) and x{{j + l)e), to = U, and tjv+i = tf. 
For a free particle the classical path is given by 

Xfree(t) = ~ X * (t - tj) + Xj (18) 

so that 



(x f ;t f \x i ;t i )=^m I H^e^^o^+i-^) 2 . (19) 

3=1 

If the particle is moving in a potential we do not know the classical path in general, however this poses no problem 
as we shall see. 

For a particle in a harmonic potential 

C = -mx 2 — —muj 2 x 2 (20) 

the classical path is 

xr — Xi cos we . 

x s ho = Xi cos uit — tij H — : smcj(£ — ti) 

sinwe 

~ x free + 0[{t-U) 2 ] 

= x free + 0{e 2 ) (21) 
We see that the classical path differs little from the free path after short durations. In general x c \ = Xf ree + 0(e) and 

V(x c i)dt= / V(x f ree + 0{e))dt^eV(x l ) + 0(e 3/2 ) (22) 



upon use of Eq. [THJ The order e 3 / 2 term does not contribute to time evolution of the system so that it is valid to 
take x c i = Xf ree between the points Xj and Xj+i and to evaluate V at some point in the interval (since the difference 
is higher order in e). The end result is that it is valid to discretise as follows 



—mi 2 — V(x) 



S c i{xj -> x j+ i) - 
The temporal evolution of a wavefunction is given by 



dt = g(x J+1 - x 3 ) 2 eV(^^) + 0(e^ 2 ) (23) 



Thus for a free particle 



,.(„/_,) = / ^ e is cl (v^) ^ {y . t y (24) 



<•(.«■:/+£)= / ^e^^) 2 iP(y;t). (25) 



The integral is washed out for large values of x — y since the exponential oscillates rapidly in that region. In fact the 
only appreciable contribution is when x — y ~ y/e. Now let 7/ = x — y and expand in e = ij 2 , rf /e, ... to get 
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ip(x;t) 



e 2 ^ 



2m ^ ' ^ 2rn 



2-Kieh 



m 



2J1 



-K 2 e 
Am 2 



■1> 



IV 



etc. 



(26) 



Thus the zero order term fixes the normalization, the first order term yields the Schrodinger equation, and the higher 
order terms are simply multiple applications of the Schrodinger equation, and so represent no new information. This 
technique can be used to obtain the quantum mechanical Hamiltonian associated with a given system. 

C. Path Integral — Quantum Field Theory. Functional Calculus 

The generalisation to field theory proceeds by discretising spacetime, ip(x) — > ipi, and defining the sum over paths 



as 



J Dp^ J Y[dipi. 



(27) 



Continuing with the development will very quickly lead to quantities like 5(ipi — ipj), d/dipi, and functions f({<Pi}) >-> 
C. The continuum version of a function like / maps elements of a function space to complex space; these are called 
a Junctionals and are denoted f[p\. A simple example would be 



(28) 



(this is the arc length functional). Manipulating functionals is called functional analysis or functional calculus. 
Unfortunately, functional methods are not nearly as well developed as real analysis or calculus, even defining functional 
integrals is fraught with mathematical complexities. But we are physicists and will happily ignore anything that gets 
in the way of progress. The functional derivative is introduced in analogy to the integral: 



i) 



ip(y) = 5(x-y). 



Now consider the problem of minimising a functional Flip} 1 , 

F[<P] = JfivM- 

Set ip = p + Sip, substitute and extract the leading change in F: 

SF = J [f(p + Sp,d^p + 5p))-f(p,d^p)} 
df 



— Sip 
op 

>1L d _ 

dp M d{d^p) 



df 



Sp. 



(29) 



(30) 



(31) 



1 Of course this is close to a physicist's heart: all classical physics follows from the minimisation of an action functional 
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The variation is arbitrary so the integrand must be zero at the minimum and we derive the Euler-Lagrange equations 
of motion (if / is the Lagrangian density and F is the action): 



df _ „ df 

Now consider the same problem from the functional point of view: upon discretising the minimum (extremum) of 
F[(f] is determined by the set of equations 



= 0. (33) 

In the continuum this is 

J^(x) m = sfifi J dy fMy)> d Mv)) = o (34) 

Taking the derivative, integrating by parts, and using Eq. [29] gives the Euler-Lagrange equations of motion. This is 
nice - the functional derivative provides a compact and transparent way to obtain the Euler-Lagrange equations of 
motion. 

For practice consider the equations of motion of (/5 4 theory: 



8tp(y) 



S[<p] 



My) 



d d x 



d d x 
1 S(p(x) 



2- 



An 



2 5<p(y) (~ d2 ^ x ^ ~ m o<P( x ) S ( x - V) - 24 ' 4 ' ^( x ) S i x ~ v) 



X 



(35) 



We now return to the path integral over fields. One can write the amplitude for a state ipi at —T to evolve to iff 
at T as 



(^(x^le-^l^-T)) = / D^%^l x) e^-^ (36) 

If one considers the limit T — > oo and assumes adiabatic switching 2 then the expression becomes manifestly covariant 
and represents the vacuum-to-vacuum transition amplitude. In the presence of an external source field, J, this is 
called the generating functional. For a scalar field theory 



Z[J] = J Dtpc iS+iJd l xJ{x) v {x) ( 37 ) 

Z is called the generating functional because it contains all of the information necessary to evaluate observables in 
scalar field theory. For example, it can be shown that the two point function defined by 



(0|T[e- l /^]|0) V ' 

is also given by 

(<Pi<P*) = ^J tpi ^ 2 e* s+i f J * | J=0 . (39) 



2 This means that the initial and final states evolve to ground states of the noninteracting system. 
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We check this explicitly at zeroth order in perturbation theory: 



f d 4 k e - ik -( x i- x 2) 

Comparing to the functional method requires knowledge of the quadratic integral 

Dlpe iff-PiM 12V2 +zfJ V = ( det M)- 1 /2 e -i//Ji(Af- 1 ) 12 J^ ( 41 ) 

Use the last equation to confirm the expression for the free propagator: 

(^W2> = ^- J T j Z [J]\j=o 

Zq 10 Ji 10J2 
= ^ ( ~Q2 i 2) 12 

We have defined the free generating functional Zq, which is obtained from Eq. [37] upon setting the interaction to zero. 
Notice that the Feynman prescription was introduced in the last step of Eq. [42j It is of course required to select the 
appropriate poles in the Feynman propagator. This can be accomplished by simply setting 

d 2 + to 2 -> d 2 + m 2 - ie 

in the action. This has the effect of multiplying the action by e~ eVT , which acts as a convergence factor on what is 
otherwise a distressingly oscillatory integrand. 

For later we note that the propagator is a Greens function of the Klein-Gordon equation: 

(dl+m 2 )A{x-y) = -iS 4 (x-y). (43) 



D. Diagrammatics 



Recall that the Gell-Mann-Low equation (Eq. [T2|) is normalized with a vacuum matrix element. This has the effect 
of eliminating all bubble diagrams from scattering matrix elements. This is desirable since bubble diagrams do not 
contribute to scattering and they introduce factors of the spacetime volume (ie, infinity) into the formalism. Diagrams 
that contain portions that do not connect to scattering particles (external legs) are called disconnected (see Fig. [!}. 




FIG. 1: (left) a disconnected diagram, (right) a connected diagram. 



In our case we must divide all expressions by Z (as we have been doing above), or equivalently, consider derivatives 
of log Z. This is most easily achieved by defining the connected generating functional: 



e ii? M =Z[J}. 

Derivatives of F automatically yield connected Greens functions. 



(44) 



FIG. 2: (left) not 1PI, (right) 1PI 



= -o— + — o— o— * - 

FIG. 3: Full two-point function in terms of the two-point 1PI diagram. 



A one-particle-irreducible (1PI) diagram (also sometimes called a proper diagram) is a connected diagram that 
remains connected when any internal line is removed (Fig. [2]). 

1PI diagrams are useful because they can be summed in simple ways to generate full diagrams. For example if an 
open circle represents a 1PI vertex and a full circle represents the analogous full diagram then the two point function 
is related to the two point 1PI diagram by a geometric sum (see Fig. [3]): 

Similarly, the full three point function is given in terms of 1PI three point functions (see Fig. 2]). 




FIG. 4: Full three-point function in terms of the three-point 1PI diagram. 



E. The Legendre Transformation 



Unfortunately, obtaining 1PI diagrams from the generating functional is not as simple as performing a normalisation. 
However the analogy to statistical mechanics helps again, and it is possible to construct a 1PI generating functional 
in analogy to the Gibbs free energy. To illustrate, consider a spin system coupled to an external magnetic field, B. 
Call the induced magnetisation M. The Gibbs free energy is defined in terms of the free energy by 



G = F~ MB (45) 

The statistical mechanical generating function is 

Z[B]=e -^ =tve -W{S)-PB-S (46) 

from which one obtains 

Now compute the change in the Gibbs free energy with respect to the induced magnetisation: 

— — B-M— 

dF dB dB 
dBdM + dM 
= B (48) 

where we have used Eq. 1471 The preferred thermodynamic configuration is determined by the minimum of G(M) 
and all quantities can be derived in terms of the induced, rather than the applied, magnetic field. 

The procedure is called the Legendre transformation in quantum field theory. In analogy, we shall define the classical 
field as the vacuum expectation value of the field operator in the presence of a source: 
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Notice that <p c i is an implicit function of the source. Computing with the free action we see 

^\x)=i[d d yA(x-y)J(y). (50) 

Using Eq. 03] we obtain 

(d 2 x +m 2 )^\x)=J(x) 

This equation tells us that it is possible to replace references of the source with references to the classical field (by 
appropriately inverting d 2 +rn 2 ). It is convenient to make this replacement implicitly via the Legendre transformation: 

r (o) b (o)j = F (o) {J] _ f ddx J{x)v (o) {x) (51) 



verify that I^ ' is independent of J. 

Ans: take jj and use the definition of ip c i. 



evaluate r(°). 

Ans: substitute , integrate by parts, and use Eq. 



You should find that the zeroth order expression for T is 

T (0) M = \ I d d x (d^ cl d»cp cl - m 2 <p 2 cl ) (52) 



Thus r^ ) is the free classical action, and we call V the effective action. 
All of this is generalised to the interacting case in a natural way: 



<p a = jjF[J\ (53) 
I> d ] = F[J]- J d d xJ(x)vd(x) (54) 
1)1 -J{x) (55) 



The last equation is obtained by taking the derivative of the second and recalling that J and tp c i are independent 
variables. As alluded to above, the effective action is a generating functional for 1PI diagrams. Specifically, the 
expansion of the effective action is in terms of 1PI n-point functions: 

°° 1 f 

Flfci] = X! ~ d d x 1 ...d d x n r( n \x 1 ,...,x n )(p c i(x 1 )---tp c i(x n ). (56) 
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F. Evaluating and Interpreting the Effective Action 

1. The Effective Potential 

If the Hamiltonian is T+V and the action is T— V, is there an analogue to the interaction, V in the effective action? 
The answer is yes, it is called the effective potential and it can be obtained by performing a derivative expansion of 
the effective action: 

rOci) = J d d X y-V ef f(cp c l) + ^Zeffiip^d^cld^ipcl + . (57) 

Notice that V e ff and Z e ff are functions of ip c i, not functionals. Thus 

V ef f (v?) = V + V 2 ip 2 + V 4V > 4 + ... (58) 

with a similar equation for Z e ff. This expression implies that it is sufficient to take ip c i equal to a constant to obtain 
the effective potential. This simplifies the 1PI expansion of the effective action (Eq. [56)) 

I> c , = const] = - / dd xi ■ ■ ■ d d x n r%,..., x n ) <{% (59) 

n=2 U ' J 

or 

T{^ cl = const) = V -f W (0, . . . , 0) ^3 (60) 
z — ' n! 

n=2 

where momentum space 1PI vertex functions have been introduced 

T(x 1 ,...,x n ) = J ^...^e^^----^^(2^) 2 % 1 + ...+p„)f( Pl ,..., Pn ). (61) 

We conclude that the effective potential is determined by Feynman graphs evaluated at zero external momentum. 

We obtain an interpretation of the effective potential by recalling that the generating functional is also a vacuum- 
to- vacuum transition element. The external source distorts all eigenstates of the Hamiltonian so we shall label states 
|n(J)) in the following. Note that if the source is time-dependent then the incoming vacuum need not be the same as 
the outgoing. We eliminate this possibility by considering J = J(x). Thus 

z[j\ = ^{M^- mJ)T \^m 
= ^E e " lB " (J)T ^/K J )> 

L J n 

( W |0(J))(0(J)|^) 
(<Pf\0)(0\<Pi) 

= e tF . (62) 

A complete set of states has been inserted in the second step and the large time limit has been taken in the third 
step. Thus (take Eq — and neglect the normalisation term) 



H - I Jipd 



F[J) = -(0(J)| 

Recalling the definition of the effective action (Eq. |55|) then gives 



|0( J)) • T. (63) 



T[y cl ) = -{Q{J)\H\Q{J))-T. 



(64) 



12 



Now let the classical field be a constant and refer to the derivative expansion of Eq. [57] to obtain 

Veff(tpci = const) = (0; (p c i\U\0; (fid) (65) 

where the notation for the states indicates that the vacuum configuration depends on the source, and hence on ip c i. 
The spacetime integral in equation[57]yields a factor of VT; dividing by the spatial volume then gives the Hamiltonian 
density that appears on the right hand side. 

Thus there is a simple and natural relationship between the effective potential and the Hamiltonian of the system. 
From this we can derive that the effective potential must be real since the Hamiltonian is Hermitian 3 . Furthermore, 
in analogy with statistical mechanics, it can be shown that the effective potential must be convex. Lastly, recall that 
ST/Sipd = — J so that the value of the (constant) classical field is determined by the simple equation 

^U=o=0. (66) 

2. The Loop Expansion 
[This section is a simplified version of chapter 11 of Ref. [l(] 

At tree order the effective potential is given by V e ff = \m^ l tp 1 cl + V((p c i). Quantum corrections can be computed in 
a clear way with the functional formalism that is called the loop expansion. 

For example, consider the computation of Z[J] for ip 4 theory to one-loop order. Let 



ip = ip c i + t) (67) 



where the classical vacuum satisfies 



ss 



ni = 0. (68) 



Expand the action 



Thus 



a s,, ,,) : / d d xC{ Vcl ) + \J v(x) 5 ,M + ■■■ (69) 



Z[J]sie i J' £ W +i J' Jw . Drjeilv^v^ ( 70 ) 



Z[J] w e 4 / c ^)+i I Jfd . det ( AJL\ r n \ 

\o<P°<Pj 

The second factor represents quantum (one-loop) corrections to the classical effective potential. 
Continuing, 

iF^i J d d x£{tp cl )+i J Jtp cl - i log det (4^) (72) 



3 Perhaps things are not so simple. Recall that scattering states are non-normalisable and hence make the Hamiltonian non-hermitian. 
Could a similar thing happen in field theory? One could, for example, interpret such an occurrence as an indication of a metastable 
vacuum. 
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and iT is the same expression without the source term. Now use the expression 

log detO = tr log O. 



(73) 



Show this. 

Ans. Diagonalise O. 



to obtain 



Substitute 



and evaluate 



r 



i , / s 2 c \ 

£(^,) + -trlog(-— j 



S 2 C 



= -0 2 -m 2 - V"(tp cl ) 



tr log(3 2 +m 2 + V"(ip cl )) 





r d d k 


VT ■ 


' (2ir) d 


VT ■ i 


[ d d k E 




J (2n) d 


VT 









log(~k 2 - 

iog(4 



+ - 4 (m 2 + V") 2 lo i 



ml + V"{y cl )) 
- < + V"(cp d )) 



V" 



A 2 



V 



(74) 



(75) 



(76) 



In the second line we have made the Wick rotation k$ — > ikd- We have taken d = 4 and regulated with a momentum 
cutoff to obtain the third line (for details on regulating field theory see section IV El) 4 . 

Remember that tp c i can be taken to be constant to obtain the effective potential. Substituting (Eq. |57|) and 
cancelling the spacetime factors gives 



V e ff((f c l) = V(ip c i) + 



16tt 2 



above 



(77) 



Unfortunately, the expression in brackets is ridiculously divergent. Dealing with divergences is called renormali- 
sation. This topic will be discussed in section IV El For now we note that the A 4 term can be neglected because it 
contributes to the vacuum energy density, which we set to zero. Also recall that mo and Ao are bare parameters, and 
must be determined by comparing predictions to experiment. Doing so permits one to replace the bare parameters in 
terms of physical parameters, called m and A. These parameters are determined by the scale at which the experiment 
is made, which we call M. The net result is 



Veff&cl) 



„-2 2 



24 



64tt 2 



{V"Y log( 



25 ,2 4 1 ,2 4 , f 2m 2 \ 



(78) 



where now V" refers to the renormalised potential, m <p + Xip 4 /24. 



4 The traditional method of computing has been presented. The observant student will note that the logarithm has units, which is 
unacceptable. Units can be supplied by subtracting logA 2 while cutting off at A. This yields a term ^A 4 log(f + (m^ + V")/A 2 ) — ► 

iA^K+V"). 
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There are some subtleties to consider here 

(i) the effective potential depends on the scale M. This scale arises when infinities are removed via the renormal- 
isation procedure. The M-dependence is absorbed into the definitions of the couplings, specifically if you choose a 
different scale, the couplings will change in such a way as to leave the potential invariant. This phenomenon goes by 
the name of the renormalisation group. 

Verify the renormalisation group. 

Ans: Show that if a different scale M' is used, the new coupling A' = A + log ^p- yields the same 
form for V e ff. 



(ii) as m — > the effective potential simplifies to 



Veff(fcl) = 



32tt 2 



&{ M 2 ' 6 



(79) 



This has a minimum at 



^«M 2 e-^. (80) 

However, the minimum is obtained by equating leading and next order (quantum) terms and thus implies that 
perturbation theory is invalid, and hence the new minimum is likely invalid 5 . 

Show that the minimum of massless ip 4 theory in three dimensions is perturbatively stable at one loop 
order. 



3. Symmetry Breaking 

Symmetry breaking happens when the ground state of a system does not have the full symmetry of the action. 
Unfortunately, the terminology is something of a malapropism as the symmetry is not really broken, rather it is 
hidden. Note also that some authors distinguish spontaneouslg broken symmetry (a scalar field acquires a vacuum 
expectation value) from dynamically broken symmetry (in which symmetry is hidden even when there are no scalar 
fields). There is a major practical difference between these: spontaneous symmetry breaking can implemented by 
simply assuming it (this is what is done in the Standard Model), while dynamical symmetry breaking requires a 
nonperturbative understanding of the system. 

Whether symmetry is broken or not, the effective potential must maintain the symmetry of the action. For example, 
in ip* theory the effective action must be even. The reason is that if it were not true, then infinities in the coefficients 
of the new terms would require renormalisation via a Lagrangian parameter that does not exist. Thus renormalisable 
field theories have symmetric effective actions. 



5 Sometimes the situation can be improved with the aid of the renormalisation group equations (see section [V El for more information). 
For the effective potential this is 

V eff (tp c i;M,X) = 0. 



», 9 a 9 9 

M h B 703 ci 

dM ' dX d<p c , 



Solving this equation and matching it to the perturbative result (recall that first order partial differential equations only give constraints 
on the solutions) gives 

f 



V. 



eff 



A 4 

24 



^Alog^/M"' 



which is called the renormalisation group improved effective action. 
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So how does symmetry breaking occur? It is possible for an even function to have two degenerate minima at, say, 
ip c i = ±v. In quantum mechanics any initial state centred on one minimum can tunnel to the other minimum and 
eigenstates have even or odd parity. Quantum field theoretic (or classical systems in the bulk limit) have infinitely 
many degrees of freedom, which one can think of as taking infinitely long to tunnel from one local minimum to 
another. Thus it is possible for an eigenstate to have mixed symmetry. This is what happens when a spin system 
acquires an induced magnetisation. So it is possible that the classical field stays in the vicinity of, say, ip c i = +v, 
giving rise to spontaneous symmetry breaking. 

The computation of the previous section begs the question: is it possible for quantum effects to break a classical 
symmetry? Although we cannot say if it happened for ip 4 theory, it is possible in general. A classic example is broken 
scalar QED, called the Coleman- Weinberg model. Ideas like this have also been used to compute bounds on the Higgs 
mass. 

Since most techniques develop the effective potential as a power series in the classical field, the simplest way to 
ensure reflection symmetry is to expand about the symmetric point ip c i = 0. Of course one could also expand about 
another point, say ip c i = v. In this case odd terms would be generated, for example T^ 3 ^ = —i\ov at tree order. But 
the sum of all terms will generate an even effective action and renormalising the even terms will lead to complete 
renormalisation of the odd terms. 

G. Goldstone's Theorem 

[This section follows the presentation in chapter 5 of Ref. Q-] 

We illustrate the power of the functional formalism by sketching a simple derivation of Goldstone's theorem. The 
theorem asserts that when a classical continuous symmetry is broken a massless particle necessarily exists in the 
spectrum. 

Consider the 0(2) linear sigma model. This Lagrangian has an 0(2) symmetry under rotations of the field: 

</>' = K(j) (81) 

where 1Z is a 2x2 rotation matrix and 



(82) 



For small rotation angles, e 



^ = l + efj 1 ). (83) 

Since the Lagrangian is invariant under this transformation (and it is on fields only), the action, the generating 
functional, and the effective action are also invariant. Thus 



Z[J] = J D<jye l 



Dcj)' q 1 S c W\+ l S J<t>' 
= J D<f>e ij ' c ^ +l i ' Jn ^ 

ps Z[J]+e I Dcj) [ (J a TT-J x a)e i f £ ^ +i ^ J ^. (84) 



The term proportional to e must vanish and we get 



— J ° - I = °' (85) 
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The log is there because we want to normalise. Now make the Legendre transformation (recall Ji — —ST/Sip^): 

Kdj a cl —)=0. (86) 

da c i Gird J 

This is a Ward identity for the sigma model that is true because of its 0(2) symmetry, regardless of whether the 
symmetry is hidden or not. 

To get the Goldstone theorem take a derivative with respect to ir c i{y) to get 

/ jap s 2 r \ 

d d x Ti c i (x) - — rT j — ^ +S(x- y) - — p-r - a cl (x) - — — — — =0. (87) 
V dn c i{y)da c i(x) oa c i{x) on c i {y)oTT c i(x) ) 

Finally, assume a hidden symmetry: 

' 7T cl \ = / 
<Tcl J \V 

One says that the symmetry is hidden (or broken) because this matrix element does not obey the 0(2) symmetry of 
the Lagrangian. 

Substituting into the last expression gives the desired result 

^ F f d d xv — = (89) 

Scr c i(y) J 8TT c i{y)5-K c i{x) 

v I d d x PT = J a (y). (90) 
J dTr c i(y)dTr c i(x) 

As the a source goes to zero one must either have that v — > (but this is the symmetric vacuum, which we do not 
consider) or 

limr„(p) = 0. (91) 

p->0 



or 



We have performed a Fourier transform and adopted a compact notation. We will show (see Eq. 11321 below) that 
the second derivative of T is the exact inverse propagator; this implies that the pion propagator has a pole at p 2 = 0, 
which means it is massless. 



[Amit, 5.11] Derive the relation 

T^{p) - T aa (jp) = -vT(p, 0, -p). 
Ans: Start with the Ward identity and take derivatives with respect to 7r c ; and a c 



H. Transition Elements — Quantum Mechanics 

The path integral is an infinite version of the familiar Riemann integral and hence should admit familiar integral 
theorems such as integration by parts. This can provide powerful nonperturbative information on quantum field 
theories. Some of this information is explored in this section. 

We start by reviewing some freshman calculus. Consider the integral 



dx e 



(92) 
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and let x = y + Af(y). One obtains 

/oo 
dy (1 + A/') e -hKv+M) 2 = j (93) 
-OO 

but /'(A = 0) = I so that the non-zero order terms in A must all vanish identically. Thus 

dye-^(J'-kyf)=0 (94) 

for any function /. This, at first sight remarkable, equation is simply a restatement of integration by parts and so is 
true for all / that diverge more slowly than 9 . 
Let's consider a more general example. 

Proceeding as before and using J dx f(x) e ax = f{4~) j dx e ax gives 

[/'(£) - + + V'AfA]l(a; k, V) = 0. (96) 

da da da da da da 

For example if V = — -|x 4 and / = 1 then 

AJ'"(a) + fcl'(a) - aJ(o) = 0, (97) 

which may be solved perturbatively in A for the behaviour of 1(a). Of course this equation also lends itself to 
nonperturbative studies. 

Let's apply these ideas to quantum mechanics. Start with the detailed expression for the path integral of Eq. [23] 

(*j\tfW t u) = /n^ e * Ef = o?(353 ' +i ^ )2 " eM=i±l±fi) ( gs ) 

and make the transformation Xj y.j + 5fj, with fj = f(jjj): 



3=1 

where x = x it x N+ i =Xf, y = yt + A/j = x it and y N+ i = y f + Af f = x f . But 

N N N 

- V 3 + A (fj+i - fj)? = ^2(Vj+i ~ Vif + 2A^(2/,+i - yj )(f j+ x - fj) + 

j=0 j=0 j=0 

2A(y l f l ~y f f f ) + 0(A 2 ,e 2 ) (100) 
where now yo = yi and i/at+i = yf so that the path integral is now in its usual form. The order A term in Eq. [Ml is 

/N . ~tf 

Dye is cl {y,y)^ fl . + ?_ {myf-V'(y)f)dt+ t ^-(y i f i -yfff)]=0 (101) 
3=1 J u 

integrating yf by parts then gives 
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N — ' r*/ 

( P >= t</, m ,m> * ) m) 

This is a remarkable and powerful equation. For example if fj = djk/e (f(y(t)) = S(t — k)) then 

(mx) = -{V), (103) 



which is the quantum analogue of Newton's equation. 

Finally consider fj = yjSjk/t or f(y(t)) = y(t)5(t — k) then 



26 

(1) = j(myy(k) + V(y(k))y(k)) 



Or 



i, (Vk+i - 2y fc + y k -i), 

>» ) 

n e 

i, /Vk+i-yk yk-yk-i\ x 

-(my k ( )) 

Um yk+1 ~ yk (y k -y k+1 )). (104) 
n e 



( (yfc+i ; yfc) (105) 

e z it 

The right hand side goes to infinity as e goes to zero so that the important paths are not those with a well-defined 
slope. In other words, infinitely spiky paths dominate the sum over paths, in keeping with the statement made in Eq. 

us 

I. Transition Elements — Quantum Field Theory 

The field theoretic analogue of the previous expressions is developed in the same way. Here we choose to work with 
the continuum notation. Consider the change of variables tp — > \ + e-F(x) m t ne generating functional of a scalar field 
theory: 

Z[J] = [ D<pe iS+i f J,fi . (106) 

Z[J] = f D X \^\ e iS\x+eF]+ifJx+iefJF 

J ox 

= J D X \5 d (x-y)(l + e^)\e iS M +i I J x (l + ie J ^-F + i e J JF 

= M i+e /f)^ l+ ' /J ^W(f F +^ 



(107) 
(108) 



The term in the first brackets above is the measure; in discretised form it is 



n^(i+^/)^n^+ e n^l; A ^ i;i ^^ x+ ^ xe<5(0) (io9) 

3 

where 1/A has been written as S(0). 
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As before, the order e term must evaluate to zero, hence 



d d x 



Sd (°) F ' ( ~T7T~\ ) + * ( lim 7^ ( "TTTT ) + J ( x ) ) F 



i8J(x) ) \y^x Sip \iSJ(jj) J ) \iSJ(x) 



Z[J] = 0. (110) 



The limit is introduced to remove ambiguity in taking derivatives. 
If F = 1 one obtains the quantum equations of motion 



((d 2 + m l) x + v'(x)) = (J)- (in) 

We will shortly see that this equation, written in terms of the effective action, is equivalent to the Schwinger-Dyson 
equations. 




J. The Schrodinger Functional 

Before moving on to a description of the Schwinger-Dyson equations, we consider one last application of the 
functional formalism. 

The path integral maintains a prominent place in quantum field theory because it is a compact representation of the 
quantisation process and because it is manifestly Lorentz invariant. Nevertheless, one can ask where, if anywhere, the 
Hamiltonian fits into a functional formalism. In fact all of the familiar language of Hamiltonian quantum mechanics 
maps over (with some modification) to quantum field theory. Thus one speaks of a Hamiltonian density, the vacuum 
wavefunctional, eigenstates, variational principles, and so on. 

The Hamiltonian for a given theory can be derived from Eq. [24] since temporal evolution is controlled by H. 
Extending to the quantum field case yields 

■9[tp;t + e) = J Drie iS ^f^[mt) (112) 

We stress that 'J is a wavefunctional that depends on a field, ip is the final field defined at time t + e; rj is defined at 
time t. Expanding the action in powers of e and recalling that 6x 2 ~ e gives 

t— *M = / d^xH^M (113) 

with 

n = + ■ v ^ + m °^ 2) + Y^- (114) 

It is tempting to interpret (ip\^) as the amplitude for the field to realise the specific value ip, but we should remember 
that only S-matrix elements carry a physical interpretation in quantum field theory. 

Let's push the analogue further and consider a variational estimate of the vacuum functional. We start with a guess 
for the free action: 

(^|* )oce'^^^ fcMfc)v(fc) . (115) 

Evaluate 

trial[ ] ~ ( 6) 
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and recall that 



(*|»7>= J Dip ty*[ip]r][[p]. (117) 
It is convenient to rewrite Hq in momentum space (and we consider d = 4): 



1 f d 3 k f 5 



K-2 



Thus 



jfDvfj^ K(2vr)^ 3 (0) + ( g 2 + ml Jfrfi) e~ f 

Evaluating the functional integral gives 

2W«] = ^(0)(27r)» J ^ (io q + ■ (120) 

Take the functional derivative with respect to the unknown lj and solve the equation to obtain 



uj k = ^k 2 + m 2 Q . (121) 

Thus the ground state functional for free scalar field theory is given exactly by a Gaussian functional with a kernel 
specified by the free particle dispersion relation. Finally, substituting gives the familiar ground state energy density: 



Use the same Ansatz and obtain the equation for the ground state wavefunctional for ip 4 theory. 

d 3 q 1 
(2tt)3 u q 



Ans: ul = k* + mt + ^S£2-± 



One final remark: Feynman has commented that the variational principle must be useless in quantum field theory 
because the set of functions probed in any variational Ansatz has measure zero. The problem above indicates that this 
is too pessimistic: we shall see shortly that the equation for oj is precisely the same as a truncated Schwinger-Dyson 
equation for the propagator. 

V. SCHWINGER-DYSON EQUATIONS 

A. The Schwinger-Dyson Master Equation 

In the functional approach, the perturbative evaluation of any n-point function for ip 4 theory proceeds from the 
expression 



Z[J]=c lV (^j) Z [J] (123) 



with V = f ■ Since Z$ is a known quadratic form, Z[J], and hence all matrix elements of products of fields, is 
computable simply by taking repeated derivatives. 

Schwinger-Dyson equations are one way (there are many) to organise the diagrams that contribute to n-point 
functions in a quantum field theory. They are popular because they naturally sum infinitely many diagrams and 
therefore automatically contain nonperturbative information. 
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The starting point is the simple identity (set F = 1 in Eq. II 10[) : 



Dip 



Sip 



= 0. 



(124) 



Extracting the functional derivative using a standard trick (see Eq. 196 j) and substituting the definition of the connected 
generated functional (Eq. |4"4"]) gives 



-iF 



SS 

5ip{x) 



5 
iSJ 



-J(x) 



(125) 



While workable, this equation is problematic because we prefer to work with irreducible n-point functions, which means 
eliminating all references to F. In general this is quite painful, involving more and more elaborate substitutions as one 
moves up the sequence of Schwinger-Dyson equations. Fortunately there is a simple trick that eliminates references 
to F from the start: the previous equation is equivalent to 



8S 
5ip(x) 



_6_ 

iSJ 



5_F 
Jj 



■ 1 = - J{x) 



(126) 



The '1' on the left hand side indicates that one should take derivatives of all quantities to the right of a derivative 
and that one stops upon taking the derivative of unity. Thus, for example, the functional derivative appearing in the 
argument of SS/Sip only acts on SF/SJ. 



Obtain Eq [126] from Eq. IT251 

Ans: work in perturbation theory or follow the strategy employed in proving the Baker-Campbell-Hausdorff 
formula. 



Now make the Legendre transformation, using the chain rule to replace the derivative with respect to the source 
with: 



iSJ 



dd 6(p ci (z) S 



iSJ 5ipd(z) 



(127) 



then replace 



S 2 F 



6<Pd(z) _ 

5 J ~ SJSJ(z)' 



(128) 



We are almost there, we only need obtain an expression for the second derivative of F in terms of the effective 
action. Consider, therefore, 



6(x — z) 



Sipci(z) 



~ j ddy 5J(x)SJ 



SJ{x)SJ(y) Sip c i(z) 



= - d d y 



S 2 F 



5 2 T 



SJ(x)SJ(y) 5ip c i(z)6ipci(y) 



(129) 



Thus 



S 2 T 



5ipd{x)Sip c i(z) 



S 2 F 



8J(x)8J{z) 



(130) 



This is the relation we need to complete the derivation of the Schwinger-Dyson master equation. Before writing this 
we first note that setting the source equal to zero implies that the vacuum expectation value of the field is also zero 



<Pcl\,J=0 = 0. 



(131) 
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(We shall discuss a more general situation below) . Finally, Eq. 11301 implies that 

<5 2 r 



S(p c l(x)S(fi c l(z) 



iPcl=0 = iA- i (x- z). (132) 



The irreducible two point function is the inverse of the exact propagator. 

We now continue with the derivation of the Schwinger-Dyson master equation; substituting Eq. 11301 into Eq. 
along with Eqs. 11281 and l54l gives 



Sep 



Pel + i 



S 2 T V 1 S 



SifclSifdJ Sip c i(z) 



1 = f- (133) 



This is the master equation equation for generating all Schwinger-Dyson equations for a scalar field theory in terms 
of irreducible n-point functions. 

Let us make this more concrete by considering (p 4 theory. Recall that (Eq. 1331) 

g = -(a 2 + ™>-^y. (134) 

Thus the master equation becomes 

- (d 2 + ml)ip cl (x) - ^ L cl {x) + A xz -\^) • 1 = - (135) 
6 \ Olfcl{z) J oip c i(x) 



Simplify using 

<5 



5(p(z) Sip{z) \6ip5tpJ xy 

We have defined the three point function 



S ( 5 2 T* \ ' 

^ x v = i 7Z7I\ ( YZIZ ) = A xa (ir azb )A by (136) 



r s3r i 

abc 6tp(a)8tp(b)6(p(c) 



Henceforth all fields will be classical fields, and we drop the cl subscript. 



prove Eq. 11361 

Ans: consider this as a matrix equation and use g(Mir') = 0. 



The final form of the master equation for (p 4 theory is thus: 



(d 2 + ml)cp x - — Up% + 3A xx (p x + A xa A xb A xc iT abc ) = - — (138) 
6 5tp x 



B. Schwinger-Dyson Diagrammatics 

Let us consider setting the source equal to zero in the master equation for tp theory, Eq. 11381 The defining equation 
[55] then implies that the left hand side is zero. If we take (p c i\j=o — as discussed above, then this implies the last 
term on the left hand side is zero, which implies that the three point function is zero 



r a f> c — 0. 



(139) 
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This is sensible: ip 4 theory is invariant under parity reflections and there should be no odd n-point functions if one 
expands about ip c i = 0. 

The Schwinger-Dyson equation for the propagator follows by taking the derivative of the master equation with 
respect to the classical field: 



= -(d x + m )5 xy - —<p x 6 xy — — - — (A xx ip x ) — — - — (A xa A xb A xc iT a bc)- (140) 

0(p X O(Py 2 2 defy 0(fy 

Notice that repeated x indices are not summed in the above; in fact the delta functions in x and y are important for 
routing momenta. However, the letters at the beginning of the alphabet are summed over! Taking the derivatives 
with the aid of Eq. 11351 setting the source equal to zero, and recognising that three terms are in fact identical, gives 



\j=o — —{d x + ml)5 xy — ^-A xx 6 xy — ^- A x bA xc A x diT a b c A ea iT d ye — A xa A. x bA. xc iT a b cy (141) 



QLp x 0LPy 2 2 



Exercise: verify this equation to first order in perturbation theory. 



It is convenient to rewrite this equation in terms of diagrams. While it is possible to develop these in terms of 
spacetime coordinates, it is useful to work in momentum space to take advantage of the translational invariance of the 
system. First we recall that the left hand side is the inverse full propagator; similarly the first term on the right hand 
side is the inverse of the lowest (tree) order propagator given in Eq. [40j The second term contains a full propagator 
that runs from x to x and a 'stub' that runs from x to y. Similarly, the last term involves three full propagators that 
run from i to a four-point vertex, whose last leg connects to the point y. Finally, we can write Eq. 11411 as shown in 
Fig.[3 



-1 


-1 


o = 





o - <t> - E) 



FIG. 5: The propagator Schwinger-Dyson Equation 



Following tradition, full propagators and vertices are represented with dots. 
We note the following: 

(i) diagram topology is specified by following indices in the Schwinger-Dyson equation 

(ii) factors of i, 1/2, -1, etc are absorbed into the definition of the diagrams except as indicated in Fig. [5j This is 
because the perturbative Feynman rules are sufficient to determine all of these factors uniquely. The explicit minus 
signs then make everything work out 

(iii) the master equation implies that there must be exactly one bare vertex in every Schwinger-Dyson equation 

(iv) the master equation restricts the form of possible diagrams that appear in Schwinger-Dyson equations 



show that 



cannot occur in any Schwinger-Dyson equation. 



Ans: use the diagrammatic language below, the starting diagrams, and the fact that there is no 'joining' 
operation. 



Although we have considerably streamlined the process of generating Schwinger-Dyson equations, taking derivatives 
gets tedious rapidly. Fortunately, the simplicity of the formalism permits an entirely diagrammatic approach. First 

we introduce a classical field as ip c i = ® . The master equation can then be written diagrammatically. One 

functional derivative yields the full propagator with source terms present, illustrated in Fig. O The notation has 
been simplified further by noting that all vertices are full except one, we therefore choose to label the bare four point 
vertex with a square. All internal lines are also full propagators. Remember that the external lines are stubs (do not 
propagate) . 



FIG. 6: Full Propagator with Source Terms 



5cp(y) 



= S(x,y) 
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8tp(y) / /~ i 



FIG. 7: Rules for Functional Derivatives. 



One can continue taking derivatives and translating the results into diagrammatic form, but it is easier and faster 
to work directly in the diagrammatic approach. The required elements are shown in Fig. [Jj 

Applying the rules to Fig. [5] yields the full equation for the three-point function given in Fig. [5] 




FIG. 8: Full Three-point Function with Source Terms. 

With the diagrammatic formalism one can proceed to quite high order. Setting the classical field to zero reduces 
the number of diagrams dramatically - we display the first three Schwinger-Dyson equations in Figs. [9l [TD1 QTJ It is 
evident that solving these equations - the first three of an infinite set - is a daunting task. Furthermore, there is no 
obvious small parameter which one can use to organise the series of Schwinger-Dyson equations. We examine these 
issues in sections IV El and IV Gl 

C. Several Fields — Alternate Schwinger-Dyson Equations 

Schwinger-Dyson equations provide a convenient methodology to sum Feynman diagrams. Of course, one need 
not sum diagrams in this way. One could, for example, sum all leading infrared divergent diagrams 6 . In fact, even 
Schwinger-Dyson equations do not admit a unique form. We illustrate this here with a simple two-field bosonic model. 

Consider two scalar fields, A and p, interacting according to the Lagrangian 



C = -d u Ad»A - -a 2 A 2 + -daipd^ip - ~mV - -gAtp 2 . 
2 2 2 M 2 2 



The field equations are 



Ja 

5S 

5p> 



-(d 2 + f i 2 )A- 



-(d 2 + m 2 )p — gAip. 



(142) 

(143) 
(144) 



Deriving the master equation requires a slight generalisation of the previous method to deal with many fields. 
Consider, therefore, a set of fields labelled ipi(x). All applications of the chain rule (etc) must sum over the index i. 
Let us change notation slightly 



<Pi(x) -> <p{i,x) -> <p{X) 



(145) 



This is what is done in the famous resolution of the infrared divergence problem in the degenerate electron gas. 
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1-7 i-7 



f- )o -<J> 



FIG. 9: The Full Propagator 



ex 



FIG. 10: The Four Point Function 



where the variable X is now understood to represent a spacetime coordinate and a field label. In this way all of the 
previous expressions remain valid, with the understanding that all coordinates are now 'generalised coordinates'. As 
a final bit of notation, if a specific field is specified then I will denote that by, eg, X tp x . 
The master equations are 



_5I_ 

5A(x) 
ST 

6ip(x) 



-{d 2 x + f i 2 )A(x) - -gtpl ~ -gA 9 „ 9x 



= -{dl + V 2 )<p{x) - gA x ip x -gA AxVx 



(146) 
(147) 



One more derivative yields the propagator equations shown in Fig. Q21 The <p field is represented as a solid line, while 
A is denoted with a wiggly line. 

There are several ways to obtain the vertex equation (which is the point of this section) . For example one can form 



S 2 T 



or 



SA Z Sip x Sip y 



_S 6 ST 

8ip x Sipy 5A Z 



(148) 



(149) 



Doing so yields distinct equations for the full vertex, illustrated in Fig. 1131 We have simplified by assuming that 
there is not mixed propagator, A l p x A y \j=o, which follows if there is no three-phion vertex. Notice that one cannot 
conclude that the diagram involving the (p* vertex is equal to that with the A 2 ip 2 vertex because the bare A(p 2 vertex 
in the other graph is in a different location in the two versions of the equation. Finally, this discussion carries over in 
precisely the same form to QED and, with extra diagrams, to QCD. 



Derive these results. 



D. Some Simple Examples 

1. ip 4 Gap Equation 

Let's rewrite the ip theory Schwinger-Dyson equation for the two-point function using Feynman rules and repre- 
sentations for the full propagator and vertex: 



A, 



F(p 2 ) 



(150) 



T 4 = T(£ u ta, t 3 ,li) (2ir) d S{li +£2 + £3+ 4). 



(151) 
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u - it- 




FIG. 11: The Six Point Function 



t 



We obtain 



FIG. 12: Phion and Aon Propagators 



F(p 2 



JT — TO, 



d d q 



2 7 (2ir) d F(q 2 ) 
d d h d d £ 2 d d £ 3 i 



,Ao 

* 6 7 (27r) d (2n) d (2n) d F{£\) F(l\) F(£j) 
i2TT) d S{h+e 2 +£ 3 -p). 



T{h,l 2 ,h,-p) 



(152) 



We see some immediate concerns: 

(i) if F ~ p 2 and T <~ const then the integrals diverge if d > 1. 

(ii) we need an expression for T. Approximating this as the bare vertex, T = — iA , is typical. This is called the 
rainbow-ladder approximation. 

(iii) we need to evaluate some pretty nasty integrals. 

(iv) we need to solve a nonlinear integral equation. 

These add up to some pretty hefty demands on our computational and analytic tool bags! 

Let's truncate heavily by ignoring the last term. Set d = 4 and simplify further by performing the Wick rotation 
to Euclidean space, 



q -> «<? 4 



and call q 2 E = q\ + q\ + q 2 + q\. Thus 



(153) 



(154) 



Notice that the integral is independent of p hence F = —p\ + const. Cut the integral off at A and call it —a to obtain 
the consistency condition: 



a = 2n £ 



A <f E dq E 



(2tt) 4 q 2 E + m 2 + \ a/2 
Evaluating the integral for large A yields a transcendental equation for a: 



(d = 4) 



(155) 



An 



y 



32-7T 2 1 + y log y 



(156) 
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with 



For small Ao on obtains 



FIG. 13: Two Vertex Equations in Aip 2 Theory. 
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2A 2 ' 



(157) 



ra + A 



32tt 2 



(158) 



If Ao > 32ir 2 then the equation admits no real solution (it has infinitely many complex solutions), implying tachyonic 
or decaying phions. Presumably this is an indication of a truncation problem or something deeper 7 . 
In three dimensions 



to 2 = m 2 + \ a -^ +0(«/AgA) . (159) 

47T v 

Recall that the coupling has units of mass in three dimensions. 

In both cases the phion mass moves to the ultraviolet cutoff, A. This is generally an undesirable situation since one 
must then cancel two large numbers (the bare mass and the cutoff) very precisely to yield a relatively small physical 
mass. This is called a fine tuning problem. Although this example is somewhat academic, the problem is very relevant 
to current quantum field theory since the Higgs mass requires fine tuning in the Standard Model. 

Finally, we make contact with the exercise in section IIV J[ in which you were asked to make a functional 
variational estimate of the ground state energy of ip 4 theory. The result was a wavefunctional of the form 
(ip\y>) ~ exp(-l/2 / ipujip) with 

^ = fc 2 + m 2 + ^ /• A J_ (160) 

4 J (2-kY uj q 

We seek the relationship of this result to the Schwinger-Dyson formalism. Since the functional approach is in three 
spatial dimensions we return to Eq. I152[ set d = 4, neglect the second term, and perform the qo integral. To do this 
note that the integral is independent of momentum and therefore F{p 2 ) = p\ — p 2 + const. Take this into account by 
setting F p = p 2 , — Lu p . Performing the qo integral then gives 



2 2 -a _2 f d q 1 , . 

p Q -u p = Po -p -m --J (161) 

which agrees with the exercise. We see that the functional variational approach yields results that agree with the 
Schwinger-Dyson equations when truncated at the Hamiltonian level. 



7 It is believed that c^ 4 theory is trivial in four dimensions. 
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2. Fermion Contact Model 



Consider a model defined by the Hamiltonian 



H 



1 



where 



d A xip(-tf- V + m)^ + - / d i xd i yi>{y)^ T a il)V{x-y)^ T a il){x) 



V(r) = ^S(r )5(r) 



(162) 



(163) 



This model represents relativistic fermions interacting via an instantaneous contact potential. The fermion currents 
in the interaction are the temporal components of a vector current "07^ V' with the addition of a flavour (or colour) 
structure denoted by the matrix T a . The coupling A is dimensionless, thus two additional powers (denoted by A) are 
required to define the interaction. Since this model is not renormalisable we take the scale A to also be the ultraviolet 
cutoff of the theory. Finally note that the Lagrangian is invariant under chiral rotations when m = 0. 

The full propagator is given by the Schwinger-Dyson equation shown in Fig. HU (this figure does not show the 
minus signs we have previously made explicit). Notice that the interaction is not dressed in this equation, as one 
might expect since the potential is a fixed function and not a dynamical quantity. 



-1 



-1 



FIG. 14: Schwinger-Dyson equation for the propagator in the fermion contact model. 
The full fermion propagator depends on two scalar functions, hence we employ the Ansatz 

S(k) 



A(p 2 )ft- B(p 2 )' 

We truncate by neglecting the four-fermion interaction in the last diagram of Fig. [14] Substituting gives 



Aft — B ft — m 



- tr(T Q T a ) 



d 4 q 



{2ir)* l0 Ai-B 



loV(p - q) 



(164) 



(165) 



Substitute for V and evaluate the Dirac trace of both sides to obtain an equation for B and multiply by ft and take 
the trace to obtain an equation for A. The result is 



A 



P 2 - l j2°F 



A( P 2 )p 2 
B(P 2 ) 

Notice that the integral in the equation for A is zero hence 



■ A r, 
A 2 



d 4 q A{q 2 )p-q 
(2tt) 4 A 2 q 2 - B 2 

d 4 g B{q 2 ) 
(2tt) 4 A 2 q 2 - B 2 



(166) 
(167) 



A = l 

Furthermore, B must be a constant, and the Euclidean space equation for B is 



(168) 



A f A d A q E 1 

B = m+ _C F BJ {27T)Ul + B2 - 



(169) 



Although this equation is quite similar to the ip 4 case studied above, the extra factor of B in the numerator makes 
an important difference. Setting m = 0, dividing by B, and evaluating the integral yields 
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where x = B j A. The term in brackets is bounded by zero and one, hence a nontrivial solution does not exist unless 

A>^. (171) 

Below this value B = 0. 

We see that there is a critical coupling above which a fermion mass is generated. Thus chiral symmetry is sponta- 
neously broken in this model. This is a nontrivial result that is impossible to obtain in perturbation theory. 

Let's return to Eq. 11671 and perform the integral over go using the appropriate contour and Cauchy's theorem. The 
result is 

B = ^% [^fdq , B . (172) 

4^ 2 A 2 J ^/ q 2 + £2 

We have regulated the integral by cutting it off at a scale A3. Notice that this amounts to a different regularisation 
scheme from that adopted in the Euclidean approach. Does this make a difference? 
Dividing by B and doing the integral gives 

1 = ( yV^T^ + \# log f 2 ) (173) 

where y = A3/A. The term in brackets is bounded by zero and y 2 so no solution exists unless 

A > 4£- (174) 

We see that selecting y = 1/V2 yields the same critical coupling as the Euclidean case. While this is reassuring, it 
is somewhat distressing that the specific dependence of B on the coupling is not the same in the two cases: physics 
should be independent of the regulator scheme chosen. The problem here is that the model is not renormalisable and 
hence residual regulator scheme dependence exists. 



3. Fermion Contact Model - Hamiltonian Approach 

We have seen in section II V J I that it is possible to apply functional methods directly to the Hamiltonian. In 
that section it was shown that the equation for the kernel of a Gaussian wavefunctional Ansatz is equivalent to the 
Schwinger-Dyson equation for the propagator when truncated at the 'Hamiltonian' level (ie, interaction terms are 
not iterated). The same observation applies to the contact fermion model considered here. Demonstrating this in 
the functional approach require a generalisation of the method to Grassmann fields, which are fields with properties 
appropriate for anticommuting quantum quantities. Since we have not discussed this (and will not), the equivalence 
to a third approach will be shown here. 

Consider a Bogoliubov transformation which is a canonical transformation in the particle basis: 

B k = cos - b k - a ■ k sin - d' k 

D k = cos - d k + a ■ k sin ^ b\ (175) 



Show that the Bogoliubov transformation is canonical. 
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The function 4>(k) is called the Bogoliubov angle and it parameterises the structure of the vacuum (recall that particle 
creation and annihilation operators arc defined with respect to an implicit vacuum). To make progress one only need 
the field contractions 



/d 3 k 
^[A + (fc)] Q/3 e ife ^ 

/r1 3 k 
{A_(k)] Pa e- ik <*-® (176) 



with 



A± = ^(1 ± sin (t)(3 ± cos (j)d ■ k). (177) 

Here j3 = 7 and a 1 = 7°7 l . 

With these rules one can compute the expectation value of the Hamiltonian with respect to the unknown vacuum. 
The result is 



f A d 3 k 

(H) = -2N C J —( S ( k )m + c(k)k) 

\ C A d 3 k d 3 n / \ 

+ 2A* CfNc J J^Wfi 1 - s{kMq) C{k)c{q)k ■ *) (178) 

with s(k) — sin0(fc), c(k) — cos(f)(k). Then the vacuum energy is minimised 



— (H)=0 (179) 



Cp f 2 7 B 



giving 

B = m + X 1TTT^ I ^dq , = = (180) 

where B = ks(k)/c(k). This is precisely the second form of the gap equation obtained in the previous section. Thus 
minimising the vacuum expectation value of the Hamiltonian is equivalent to solving the Schwingcr-Dyson equations 
when truncated at rainbow-ladder order. 



4. Ladder QED 

[This section follows the description in Ref. Q.] 

As a final example, we consider quantum electrodynamics in four dimensions in rainbow ladder approximation. The 
Lagrangian is 



C = -jF^Fp, - i(<9 • A) 2 + ip(id ■ 7 - m)il> - e^A ■ ji/j. (181) 

The parameter £ is the (covariant) gauge fixing parameter. The Schwinger-Dyson equation for the electron propagator 
is then given by Fig. [15] The Schwinger-Dyson equations have been truncated by employing the bare photon 
propagator given by 

D »- = STi:(^-< 1 -f'^)- < 182 » 

The gauge parameter is often taken to be zero since in this gauge £ itself is not renormalised (ie, it stays at zero). 
This is called Landau gauge. 
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FIG. 15: Gap Equation for QED in the rainbow-ladder approximation. 

Following the same procedure as in the previous examples we can derive that A — 1 in Landau gauge, and hence 

d 4 q (3 + 0)B{q 2 ) 



B(p ) = m — ie 



(27T)4 {q 2- B 2 ){p-qf 



(183) 



(the zero above comes from £). 

Remarkably, the angular integral has a very simple form: 



(p-q) 2 



2tt 2 



l -6{q >p) + ^6{p > q) 



P 



(184) 



Substituting, setting m = 0, Wick rotating, and regulating with a cutoff yields 



B(-P%) 



3f_ 

8tt 2 



1 

Pe Jo 



■ 2 i d i<i 3 - 



B(q) 



B 2 



dqq- 



B 2 



(185) 



This is a relatively simple integral equation to solve numerically since there are no angular integrals to evaluate 
and the p^-dependence is easy to implement. However, the equation can solved analytically by converting it into a 
nonlinear differential equation. 

Taking the derivative with respect to p\ gives 



4 , 3e 2 f PE , , B 
PeB' = / dqq 



hence 



8n 2 Jo ^ q 2 + B 2 



(p%B'y = -^ P i^- 



lfa 2 ^p 2 E + B 2 



(186) 



(187) 



This nonlinear differential equation is equivalent to the integral equation once we have specified the boundary 
conditions that are implicit in the integral. To obtain them, substitute Eq 11871 into Eg 11851 and integrate. The result 
is 



(p 4 B')l P =o-0 



(188) 



and 



(B+p 2 B')\ p=A = (189) 
We proceed by noting that Eq 11871 is homogeneous for large momenta and hence has the simple power law solution 



B p -l±Vl-3e 2 /(4^) 

This satisfies the first boundary condition, but can only satisfy the second if the power is imaginary. 



(190) 
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show this. 

Ans: if the power is imaginary, B is oscillatory and hence it can be arranged that (B +p 2 B r ) = at p = A. 



Thus a nontrivial gap function is generated (in rainbow ladder approximation in Landau gauge) if 

a>a* = ~. (191) 

This is a surprising result. In a world described by quenched QED, massless electrons in the field of a heavy nucleus 
with charge Z > 144 would gain a mass. 

Before getting carried away it is worth observing that 

(i) this result was obtained in Landau gauge, one should confirm that it remains accurate in other gauges. A clue 
is obtained from Eq. 11831 the factor 3 is replaced by 3 + £ in general and naively (naive because we ignore A) a+ 
changes to 7r/(3 + £). Clearly this is not a satisfactory situation. The problem, of course, is lack of gauge invariance 
induced by our simple rainbow-ladder approximation. This can be addressed by using a more complicated vertex 
Ansatz 8 or by solving the Schwinger-Dyson equation for the electron-photon vertex. 

(ii) the mass is not given by B or B / A since these are Euclidean quantities. One must first rotate back to Minkowski 
space and then solve the equation A(p 2 )p 2 ~ B(p 2 ) = for p 2 = m 2 . 

E. Renormalisation 

Except for a few warnings, we have avoided the issue of renormalisation thus far. This is an enormous topic and 
we cannot possibly do it justice here. Nevertheless, it is important to see how renormalisation fits into the functional 
methodology and Schwinger-Dyson equations. Thus a quick review is in order. 

Renormalisation consists of three steps: 

A. regulate 

B. renormalise 

C. remove scheme dependence. 

1. Regulation Schemes 

Field theories are rendered finite by introducing a parameter. This can be done in any way, preferably in a way 
that preserves the symmetries of the action 9 Some popular regulation schemes are 

(i) dimensional regularisation. Momentum integrals are made finite by computing in a complex dimension d. 
Continuing to d = 4 then reveals infinities of the type l/(d— 4). This is the most popular regulation scheme due to 
its (relative) simplicity and symmetry preserving nature. 

(ii) Pauli- Villars. Integrals are made finite by introducing fictitious particles that contribute with opposite sign. 
This scheme also preserves symmetries but can become cumbersome. It is useful in dealing with ambiguities in the 
definition of 75. 

(iii) lattice. A spacetime discretisation is employed. If the lattice spacing is a and the lattice size is La this cuts off 
ultraviolet momentum integrals at 1/a and infrared integrals at 1 /(La). It is very useful for numerical methods. 

(iv) momentum cutoff. Simply cut off momentum integrals. This scheme violates translation invariance symmetry, 
which can sometimes be inconvenient (but is not a show stopper). 



The Ansatz is called the Ball-Chiu vertex. This helps but does not solve the problem. We note that Landau gauge can be regarded as 
preferred since A = 1 brings the Ball-Chiu vertex near to the bare vertex. 

This is not mandatory: as long as a symmetry breaking regulator does not change the universality class of the theory, the continuum 
limit can be recovered. 
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2. Renormalisation Schemes 

Recall that the bare parameters of the theory are not physical. Their values must be adjusted so that the predictions 
of the theory agree with experiment. This is true for a completely finite theory and remains true for field theories with 
divergences, with the simple extension that bare parameters can be infinite when the regulator is removed. Thus, 
infinities are absorbed into the bare parameters of the theory. One would expect that this is generally impossible 
since infinities of many types can be generated by loop diagrams, while there are only a few parameters present in any 
theory. Remarkably, the program is possible in a wide class of field theories, which are hence called renormalisable 
field theories. Nonrenormalisable theories require the introduction of an ever growing list of parameters with which 
to absorb infinities 10 . 

In practice 'absorb infinities' means define a set of renormalisation conditions that the theory must obey. For 
example, the mass of the electron should be 0.511 MeV and the value of sola / d£l(e + e~ — > 77) should be 25 GeV 2 
nb/sr at cos 0cm = 0.28 and Ecu — 29 GeV. Renormalisation conditions are chosen to match the problem at hand, 
and can sometimes be quite abstract. 

In general there must be as many renormalisation conditions as there are monomials in the theory. For ip 4 theory 
this means three conditions. These specify the couplings m , A , and the strength of the field (corresponding to the 
term l/2(dip) 2 ). Lastly, all renormalisation schemes introduce another scale into the theory, called the renormalisation 
scale or the subtraction point. 

(i) momentum subtraction (MOM). Fix the residue of the propagator to be unity. Fix the pole of the propagator to 
the physical mass, T^(p 2 = to 2 ) = 0. Fix the scattering amplitude at one point, eg, A = —i\ at s = t = u = 4/3to 2 . 
In order, these conditions fix tpo, mo, and Ao. MOM is the most transparent scheme for scattering processes involving 
physical particles 11 . 

(ii) minimal subtraction (MS). Directly absorb infinities in bare parameters. 

(iii) modified minimal subtraction (MS ). Like MS but absorb some constants as well. 

(iv) fastest apparent convergence (FAC). Choose the renormalisation scale and scheme so that the last term in the 
perturbative expansion of a physical quantity is zero. 

(v) minimal sensitivity (PMS). Choose the renormalisation scale and scheme so that the prediction for a physical 
quantity is maximally insensitive to variations in the scale. 

The appearance of the renormalisation scale can be obscure. For example, in MOM one can set the scattering 
to a measured value at a general energy, M. In dimensional regularisation the coupling must be replaced with 
Am = Ao-/l/( 4 ~ d ' to obtain the correct dimensions (see Eq. |9|). Once the renormalisation conditions have been 
imposed it is possible to remove the regulator scale A, leaving only the renormalisation scale M. 

3. The Renormalisation Group 

It is rather dispiriting to work hard and end up with an expression for a scattering cross section that depends on 
an arbitrary scale M. Naively, your result can be any numerical value and nothing has been learned. But things are 
not so bad. Recall that to and A have been determined by comparing to experiment at the scale M. Thus we imagine 
to = m(M) and A = X(M). Predictability arises because changes in M lead to commensurate changes in to and A 
such that the predictions remain invariant. This behaviour is called the renormalisation group and can be determined 
by imposing 

-^(physical quantity) = 0, (192) 

called the renormalisation group equation. We have seen the application of this method in section IIVFI 

If all of this is getting overwhelming, there is another way to think of renormalisation. Regulate a theory with a 
scale, A. Fix A and evaluate several physical quantities and adjust the parameters of the theory so that the predictions 
agree with experiment. Call these too (A) and Ao(A). Change A to A' and repeat. You will trace out the functions 
too (A) and Ao(A) which reproduce experiment. Now compute some other quantities using these functions. If the 
theory is renormalisable, they will be stable as A is varied. And if your theory is right and renormalisable it will not 
matter which observables you use to determine toq and Aq. 



In the old days this was regarded as anathema. These days it is acceptable to deal with nonrenormalisable field theories if the parameter 
list can be ordered in some way, for example by energy. Theories in this class are called effective field theories. 
A physical particle is one that can reach a detector, like an electron. A quark is not a physical particle. 
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4- Renormalising the Effective Potential 

A scheme like MOM fits reasonably well into the functional formalism since it imposes conditions on the two- and 
four-point functions. To make contact with the effective potential we perform a derivative expansion (Eq. 1571) 

rfo>) = j d d x (-V ef f(<p) + \z efJ {fp)d^ip + ...). (193) 

Thus 

V eff (ip) = V + V 2 <p 2 + V i ip 4 + ... (194) 
with a similar equation for Z e ff. With this expansion the three MOM conditions are analogous to 12 



Z eff (0) = Z = 1 



'2 



±^ Veff \ v=0 = 2V 2 = rh 2 

^K//| v =o = 4!V 4 = A. (195) 

The tildes indicate that these renormalisation conditions need not refer to m and A of MOM. Alternatively, if the 
limit m — > is desired it is useful to use 

^-j V e ff\ v =M = Am (196) 

since this avoids infrared divergences that appear in the effective potential in the massless limit. This condition was 
used to derive the renormalised effective potential given in Eq. [751 



5. Renormalising Schwinger-Dyson Equations 

The first task in renormalising Schwinger-Dyson equations is regulating. Since one must employ numerical methods 
dimensional regulation is out of the question. This leaves Pauli-Villars, which is awkward, and momentum cutoff. 
Implementing a cutoff is numerically simple, but one must remember that it violates translation invariance and 
therefore care must be taken. For example, the computation of the vacuum polarisation tensor (the photon self 
energy) with a momentum cut off, A, generates a term proportional to h 2 g^ v . This term violates gauge invariance, 
but can easily be isolated and does not cause fundamental problems. 

There are several ways to impose renormalisation conditions on the Schwinger-Dyson equations. If you are familiar 
with it, using the perturbative method of counterterms works well. Another method involves forming judicious 
subtractions. Consider, for example, the two point function of ip 4 theory given in Eq. 11521 We require that the pole 
location is at the physical mass, p 2 = m 2 , which implies 



F{p 2 



Now subtract zero from both sides of 11521 



(197) 



F{ P 2 



2 2 

p — m 



d d q 



2 J (2ir) d F(q 2 ) 
.A f d d t x d d £ 2 1 1 
If J (27r) d (27r) d F{1\) F{1\ 



T{£ 1 ,£ 2 ,p-£ 1 -e 2 ,-p) T{l u l 2 ,p-h-t 2 ,-p). 



d d q 1 
(2^f(?) 



F((h+i 2 -pr) 



F((h 



pY) 



(198) 



The exact mapping would require Fourier transforming V2, V4, and Zq. This prescription is a little simpler. 
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Notice that the 'bubble' diagram drops out of the equation. This makes sense, it only contributed a constant to F, 
which can immediately be absorbed into the definition of m. If we had truncated the Schwinger-Dyson equation at 
this point (ie, not included the second diagram) we would be done. 

Take d — 4 and consider the degree of divergence of Eq. 11981 The four point function is dimensionless and so can 
at worst diverge like log A, so we ignore it. At worst F(£ 2 ) ~ £ 2 so in the ultraviolet region the integral behaves like 



d i £ 1 d i £ 2 



1 1 

J2J2 



P- 



(p 2 - to 2 ) log A. (199) 



The degree of divergence has been knocked down from A 2 to log A. Our expression for F now looks like 

F{p 2 ) = (p 2 - to 2 )G(A/ p ) (200) 

with 

A 2 

G(A/p) ~ A 2 (log + finite). (201) 
p- 

The logarithmic divergence is eliminated by imposing the residue condition 

^F(p 2 )| p2=m2 = 1. (202) 

Thus G(A/m) must be unity. Subtract = G(A/to) — 1 again to obtain the renormalised Schwinger-Dyson equation 
for the two point function: 

F(p 2 ) = {p 2 - m 2 ) ■ [G(A/p) - G(A/to) + 1] (203) 

where G is determined by Eqs. 11981 and 12001 Notice that the dependence on the cutoff is eliminated since the 
A-dependent portion of the term in brackets is 

A 2 A 2 to 2 . 

log— -log— = log— . 204 

pZ TO p z 

A single subtraction will suffice to render the equation for finite and will replace references to Ao with A. Since 
tp theory is renormalisable in four dimensions, this is all that is required to make all other Schwinger-Dyson equations 
finite. 



F. Truncation Schemes 



Although Schwinger-Dyson equations are equivalent to the effective potential, and hence to the generating func- 
tional, and hence to the quantum field theory, and they permit resumming infinitely many Feynman diagrams, they 
still form an infinite hierarchy of equations. Thus some form of truncation is required. Truncation amounts to 
searching for a small parameter. Small parameters that have been employed in the literature are 

(i) coupling constant. If A is small diagrams with many vertices are suppressed and the Schwinger-Dyson equations 
can be ordered. Of course this is equivalent to standard perturbation theory, and one is advised to use well-established 
perturbative methods in this case. 

(ii) 1/N, l/N c . If the number of fields (either labelled by an external index or an internal degree of freedom such 
as flavour or colour) is large, diagrams that contain fewer fields are suppressed, supplying an ordering mechanism. 

(iii) classical limit. If classical physics dominates the system one can order diagrams in powers of H. This method 
generates the loop expansion of the effective action. 

(iv) high density or temperature. If a system is coupled to a heat bath, thermodynamic quantities such as temper- 
ature and density (high or low) will sometimes provide a useful ordering principle. This method is used, for example, 
in the hole-line expansion for low density nuclear matter. 
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(v) infrared dominance. Many systems are dominated by long range interactions. Models of these systems typically 
contain (naive) infrared divergences. In these cases, ordering diagrams according to their degree of infrared divergence 
can yield a useful truncation scheme. This method was employed in the resolution of the infrared divergence problem 
in Fermi gasses by Brueckner and Gell-Mann. 

(vi) guess. Guesses typically take the form of truncating the Schwinger-Dyson equations at low n where n refers 
to the n-point function. While this is an organisational scheme, it is not a truncation scheme. Unfortunately this 
method is the one most commonly employed in the literature. It is not necessarily bad, but correctly implementing 
it requires that one check the efficacy of the guess and then confirm that the guess is robust. 



G. Numerical Methods 



In general it is impossible to solve Schwinger-Dyson equations analytically and one must resort to numerical meth- 
ods. Unfortunately there is no general method for solving nonlinear coupled integral equations and one must live by 
one's wits. Strategies that have been attempted are 

A) reduce the functional problem to an algebraic problem. This is done in one of two ways 

(i) expand functions in terms of a convenient basis: A(p) = CiTi(p), where hopefully the properties of T allow 
some integrals to be performed exactly (this is only the case for very simple kernels, which means highly truncated 
Schwinger-Dyson equations). A benefit of this approach is that it is simple to evaluate quantities like A(|p + <?|), which 
are often required. A demerit is that the basis may be too restrictive to permit an accurate solution (for example, 
obtaining the oscillatory solutions of rainbow-ladder QED would be very difficult with polynomial basis functions). 

(ii) place the functions on a momentum grid: A(p) — > A4 = A(pi). This method has the benefit that integrals are 
easy to obtain (numerically) and that the solution is not confounded by basis assumptions. Alternatively, it makes 
evaluation of A on non-grid points difficult. 

B) solve the resulting algebraic equations for the Ai or c\. This is a difficult problem in general. Possible approaches 
are 

(i) iteration one could attempt this by simply iterating the gap equations, using the LHS from the previous iteration 
to evaluate the RHS. This is the simplest possible procedure but suffers from instability, as is easily verified in one 
dimension. 

(ii) iterative Newton- Raphson stability can be improved by linearising the gap equations to obtain a new estimate 
of the root. Thus if one wishes to solve 



Xi = fi{{x}) (205) 
expand about an assumed solution and neglect higher order terms 

h 



Xi + 8xi = fi({x}) + j^lajx, (206) 
axi 



and solve for the corrections 



£ ( s » + 7^7; • ) **j ■' ■ /<({*»• ( 207 ) 



j 



One then sets X4 = Xi + 8xi and iterates. In general this algorithm converges to a sensible solution more often than 
simple iteration. 

(iii) minimisation construct the function G({x}) — ^2Axi — fi({x})) 2 and minimise it. This can be a powerful 
method but suffers from the bane of all minimisation problems, that of finding the global minimum in a multidi- 
mensional sea of local minima. In fact, G often has local minima in the neighbourhood of zeroes of the functions /j. 
Unfortunately, these need not be at all close to solutions to Eq. 12051 

In fact all of the techniques discussed here will fail miserably unless one starts very close to the solution. Thus it 
is vital that the practitioner not abandon theoretical investigations too early. One must carefully track and deal with 
singularities in the equations, understand asymptotic behaviour, and develop decent analytic approximations to have 
any hope of obtaining reliable numerical solutions. And it is best to be prepared for a lot of heartbreak. 
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VI. CONCLUSIONS 
A. Summary 

I hope that this short introduction has convinced you that functional methods are a very powerful and flexible 
tool for examining quantum field theory. At the least, they are more compact than the canonical approach, although 
somewhat more 'formal'. Dealing with advances topics such as curved spacetime, anomalies, or Faddeev- Popov gauge 
fixing is also much easier in the functional approach. Lastly, we have seen how the familiar machinery of quantum 
mechanics can be carried over to quantum field theory in the Schrodinger functional approach. 

Schwinger-Dyson equations are a powerful method for organising and summing infinitely many Feynman diagrams. 
As such they represent one of the few ways to obtain nonperturbative information on a field theory. And we have 
seen that in some cases it is very easy to obtain such information. The equivalence to other nonperturbative methods, 
the Bogoliubov canonical transformation and the Gaussian variational Ansatz, has also been shown. However, it is 
clear that the Schwinger-Dyson method is much more general than either of these. It is also easy to generalise the 
formalism to finite temperature and density. 

Nevertheless, problems lurk. The technical challenge in solving coupled nonlinear integral equations over many 
variables can be daunting. More fundamentally, deciding how to truncate the equations is vital and, in most cases, 
not obvious. There is room for improvement here! 

B. Where to go from here 

These notes are an introduction to an enormous field with applications across physics; much has been omitted! For 
starters fermions were almost entirely neglected due to spacetime limitations. A large literature exists on properties 
of path integrals, since their definitions are subtle. Applications to the Standard Model and attempts to go beyond 
the Standard Model abound. Thus, I include a small bibliography to help explore these topics. 

1. Original Articles 

P.A.M. Dirac, Physikalische Zeitschrift der Sowjetunion, 3, 64 (1933). The germ of the path integral. 
R.P. Feynman, Rev. Mod. Phys., 20, 367 (1948). Feynman develops the path integral. 

R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, New York, 1965). Appar- 
ently this book has become a collector's item, as new copies are going for $600 on amazon.com! Look for the rare 
ones that forgot the word 'Path' in the title. 

FA. Berezin, The Method of Second Quantisation, (Academic Press, New York, 1966). Introduced the use of 
Grassmann variables to deal with fermions. 

S. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973). Effective potentials and hidden symmetry. 

2. Books on Functional Methods 

R.J. Rivers, Path Integral Methods in Quantum Field Theory, (Cambridge University Press, New York ,1990). 
L.S. Schulman, Techniques and Applications of Path Integration, (Dover, London, 2005). 

M.S. Swanson, Path Integrals and Quantum Processes, (Academic Press, New York, 1992). Nope. Not related. 

3. Other Books 

D.J. Amit, Field Theory, the Renormalization Group, and Critical Phenomena (World Scientific, Singapore, 1989). 
A nice introduction, with a rapid approach to advanced topics, applications tend to condensed matter physics. 

P. Ramond, Field Theory: a Modern Primer, (Bcnjamin/Cummings Publishing, Reading, MA, 1981). A very nice 
pedagogical introduction to field theory and the functional method. 

VA. Miransky, Dynamical Symmetry Breaking in Quantum Field Theories, (World Scientific, Singapore, 1993). 
Discusses dynamical symmetry breaking with a focus on chiral symmetry breaking in a number of systems, including 
superconductivity, the Nambu-Jona-Lasinio model, QCD, the Standard Model, and tchnicolour. Chapter 7 is an 
excellent introduction to the effective action. 
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C. Itzykson and J.-B. Zuber, Quantum Field Theory (McGraw-Hill, New York, 1980). The standard reference. 
Unfortunately, functional methods are scattered throughout the text. Look in sections 5-1, 6-2-2, 9-1, 10-1, and 
11-2-2. Treatment of the Schwinger- Dyson equations is light. 

V.N. Popov, Functional Integrals and Collective Excitations (Cambridge University Press, New York, 1987). Ap- 
plications to condensed matter physics; does not make heavy use of the functional approach. 

M.E. Peskin and D.V. Schroeder, An Introduction to Quantum Field Theory (Westview Press, 1995). See chapter 
11. 



4- Review Articles 

Applications to QCD and other field theories are reviewed in many articles. The following will get you nearly up 
to date. 

R. Alkofer and L. von Smekal, The infrared behavior of QCD Green's functions: Confinement, dynamical symmetry 
breaking, and hadrons as relativistic bound states, Phys. Rept. 353, 281 (2001). 

C. S. Fischer, Infrared properties of QCD from Dyson- Schwinger equations, J. Phys. G 32 (2006) R253. 

C.K. Kim and S. K. You, The Functional Schrddinger Picture Approach to Many-Particle Systems, 
|cond-mat/0212557| 

P. Maris and C. D. Roberts, Dyson-Schwinger equations: A tool for hadron physics, Int. J. Mod. Phys. E 12, 297 
(2003). 

M. Pennington, Swimming with Quarks, hep-ph/0504262 

C. D. Roberts and S. M. Schmidt, Dyson-Schwinger equations: Density, temperature and continuum strong QCD, 
Prog. Part. Nucl. Phys. 45, SI (2000). 

C. D. Roberts and A. G. Williams, Dyson-Schwinger equations and their application to hadronic physics, Prog. 
Part. Nucl. Phys. 33, 477 (1994). 
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